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Abstract 

A weak-rate library aimed at investigating the sensitivity of astrophysical environments to variations 
of electron-capture rates on medium-heavy nuclei has been developed. With this library, the sensitivity 
of the core-collapse and early post-bounce phases of core-collapse supernovae to nuclear electron-capture 
is examined by systematically and statistically varying electron-capture rates of individual nuclei. The 
rates are adjusted by factors consistent with uncertainties indicated by comparing theoretical rates 
to those deduced from charge-exchange and /3-decay measurements. To ensure a model independent 
assessment, sensitivity studies across a comprehensive set of progenitors and equations of state are 
performed. In our systematic study, we find a -|-16/-4% range in the mass of the inner-core at the 
time of shock formation and a ±20% range of peak i/e-luminosity during the deleptonization burst. 
These ranges are each five times as large as those seen from a separate progenitor study in which we 
evaluated the sensitivity of these parameters to 32 presupernova stellar models. It is also found that the 
simulations are more sensitive to a reduction in the electron-capture rates than an enhancement, and 
in particular to the reduction in the rates for neutron-rich nuclei near the W = 50 closed neutron-shell. 

As measurements for medium-heavy {A > 65) and neutron-rich nuclei are sparse, and because accurate 
theoretical models which account for nuclear structure considerations on the individual nucleus level 
are not readily available, rates for these nuclei may be overestimated. If more accurate estimates 
confirm this, results from this study indicate that significant changes to the core-collapse trajectory 
can be expected. For this reason, experimental and theoretical efforts should focus in this region of the 
nuclear chart. 


1. INTRODUCTION 

The study of interactions mediated via the weak nuclear 
force is of importance to a large number of fields in physics. 
However, it is of particular importance to the field of as¬ 
trophysics because of the longer timescale on which weak 
interactions operate as compared to the strong and elec¬ 
tromagnetic interactions. This is evidenced by the impact 
that new insights into weak reaction physics have on as¬ 
trophysical models (Langanke & Martmez-Pinedo 2003). 
Specifically, electron-capture reactions play a prominent 
role in high-density environments such as those found in 
the late stages of massive star evolution (Heger et al. 2001; 
Martmez-Pinedo et al. 2014), thermonuclear (Iwamoto 
et al. 1999; Brachwitz et al. 2000) and core-collapse su¬ 
pernovae (CCSNe) (Hix et al. 2003; Janka et al. 2007), 
neutron stars (Gupta et al. 2007; Schatz et al. 2014), and 
compact object merger events (Goriely et al. 2015). Re¬ 
alistic simulations of these environments rely on accurate 
nuclear physics inputs including electron-capture rates. 

Electron-capture rates depend sensitively on Gamow- 
Teller (GT) transition-strength distributions in the /3+ 
direction. These transition strengths characterize nuclear 
excitations in which a single unit of spin and isospin are 
transferred {AS, AT = 1), with no transfer of orbital 
angular momentum (AL = 0). While the main com¬ 
ponent of electron capture occurs on the ground state 
configuration of a nucleus, in high temperature stellar 
environments electron captures on thermally-populated 
excited states of the parent nucleus can also contribute 
significantly to the overall rate (Langanke & Martmez- 


Pinedo 2000). Unfortunately, it is difficult to obtain 
information about transitions from excited states in the 
laboratory. Gompounding the problem is the fact that in 
order to accurately include electron capture in simulations 
one must include electron captures on a wide range of 
nuclei. Hence, in general one must rely on theoretical mod¬ 
els for a complete description of stellar electron-capture 
rates. On the other hand, measurements of Gamow-Teller 
strength distributions in a representative set of nuclei 
are important for the development and benchmarking 
of robust theories. At the same time, it is critical that 
theoretical and computational efforts provide guidance to 
experimenters on which measurements to perform. 

Presently, configuration-interaction (shell-model) cal¬ 
culations are the primary method for producing reliable 
GT strength distributions near stability in the sd- and 
pf- shells (8 < [N, Z] < 20 and 20 < [N, Z] < 40, respec¬ 
tively) for electron capture on both ground and excited 
states (Oda et al. 1994; Langanke et al. 2003). Quasi¬ 
particle random-phase approximation (QRPA) calcula¬ 
tions have also been utilized to estimate GT strengths 
for large sets of nuclei, but only where transitions from 
the ground state are considered (Paar et al. 2009; Moller 
& Randrup 1990; Dzhioev et al. 2010; Niu et al. 2011; 
Nabi & Klapdor-Kleingrothaus 2004). Furthermore, com¬ 
prehensive sets of electron-capture rates (as a function 
of density and temperature) for a large number of nuclei 
based on QRPA calculations have not been published. 

Direct and indirect experiments, such as /3-decay and 
charge-exchange (CE) measurements respectively, pro- 
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vide robust benchmarks for theoretical GT strengths and 
therefore are crucial for our understanding of astrophysical 
electron-capture rates. Unfortunately, electron-capture 
and /3-decay experiments can only access states in a lim¬ 
ited Q-value window. Furthermore, for neutron-rich nuclei 
/3-decay only provides information in the j3~ direction, 
which is of limited use for electron-capture studies. In¬ 
termediate energy (> 100 MeV/u) CE reactions in the 
j3~^ direction, however, connect the same initial and final 
states as electron capture, providing information about 
transitions up to high excitation energies, and are thus 
well suited to study the full Gamow-Teller strength dis¬ 
tribution of interest. At these energies, CE measure¬ 
ments are accurate at the ~10% level and are therefore 
able to provide rigorous tests of theoretical Gamow-Teller 
strengths and derived electron-capture rates. 

Recently, the results from (n,p), (d,^He), and (t,^He) 
CE reactions on nuclei in the p/-shell were systemati¬ 
cally compared (Cole et al. 2012; Scott et al. 2014; Noji 
et al. 2014) with shell-model calculations using the KB3G 
(Poves et al. 2001) and GXPFla (Honma et al. 2005) effec¬ 
tive interactions in the p/-model space, and with calcula¬ 
tions based on the QRPA formalism of Moller & Randrup 
(1990). The authors compared shell-model and QRPA de¬ 
rived electron-capture rates against those derived from CE 
measurements. It was found that the QRPA calculations 
systematically overestimate the electron-capture rates 
(^100-3000%, depending on density and temperature), 
whereas the shell-model estimates produce rates similar 
to those measured experimentally (~l-50%) (Cole et al. 
2012). Unfortunately, shell-model calculations are compu¬ 
tationally challenging for nuclei beyond the p/-shell, and 
therefore weak rates used in high-density astrophysical 
calculations most commonly rely on less accurate meth¬ 
ods. In each of these cases, systematic and random error 
exist, and it is therefore important to understand the 
sensitivity of astrophysical simulations to uncertainties 
in these rates. 

Sensitivity studies are useful tools for guiding theoreti¬ 
cal and experimental efforts because they highlight nuclei 
that should be given particular focus, and they indicate 
the accuracy with which the parameters of interest need 
to be known. They also illustrate how strongly the cur¬ 
rent parameter uncertainties affect the outcome of the 
astrophysical simulations. In this work, we perform ~ 150 
collapse simulations examining how systematic and statis¬ 
tical variations of the electron-capture rates impact the 
collapse, bounce and pre-explosion phases of core-collapse 
supernovae simulations over a range of presupernova pro¬ 
genitors and equations of state (EOS). We describe the 
development of a modular and open-source weak reaction 
rate library for use in astrophysical simulations and its 
first implementation using the stellar core-collapse code 
GRID (O’Connor & Ott 2010). In the following sections, 
we show that the inner core of the protoneutron star (PNS) 
and the observable peak neutrino-luminosity from core 
bounce and shock formation depend sensitively (-I-16/-4 % 
and ±20%, respectively) on the electron-capture rates of 
neutron-rich nuclei. As variations on this level are not 
easily reproduced from uncertainties in other inputs to 
the simulations, they motivate the development of new 
theoretical models for electron-capture rates as well as rel¬ 
evant measurements, which together will constrain these 
and other key parameters discussed in this work. 


In Section 2 we discuss previous development of weak 
rates for astrophysics, as well as the implementations 
made in this work. In Section 3 we motivate the impor¬ 
tance of electron-capture rates during the core-collapse 
phase of CCSNe, and in Section 4 we describe the codes 
we have utilized and developed for this work. We de¬ 
tail the sensitivity studies performed in Section 5 and 
conclude in Section 6. 

2. ASTROPHYSICAL WEAK INTERACTION RATES 

From the birth to the death of massive stars, the weak 
nuclear force is a primary actor in the story of stellar 
evolution (Langanke & Martinez-Pinedo 2003). Weak 
interactions are important ingredients for nucleosynthesis 
and also for the internal structure of evolving stars, as they 
sensitively determine the electron-to-baryon ratio Ye and 
the iron-core mass just prior to core-collapse (Heger et al. 
2001). Unlike the conditions present during quasistatic 
stellar evolution, however, in the core of a collapsing star 
the density and temperature are high enough that nuclear 
and electromagnetic reactions equilibrate (Iliadis 2007). 
Weak reactions, on the other hand, operate much more 
slowly and thus continue to affect the nuclear composition, 
the neutrino emission, and ultimately the dynamics of 
the entire event. 

As compared to other semi-leptonic weak interactions, 
electron capture has a particularly remarkable impact 
on the core-collapse environment (Langanke & Martnez- 
Pinedo 2014). In the final stages of a star’s life, the 
nuclear-energy generation rate of the core that normally 
sustains a star against gravitational collapse is absent be¬ 
cause the core is composed of highly stable iron peak nu¬ 
clei. Instead, at these late times the electron-degeneracy 
pressure provides the primary stability against collapse. 
It is therefore apparent that electron captures that re¬ 
move electrons from the system will have dramatic conse¬ 
quences for this environment. Furthermore, the electron 
chemical potential is sufficiently large to overcome 
Q-value restrictions, and so the electron-capture rates are 
significant. 

Just prior to and during the early moments of collapse, 
other weak interactions can also play a role. Martinez- 
Pinedo et al. (2000) have shown that /3“-decay can tem¬ 
porarily compete with electron capture when Ye=0.42- 
0.46, which can occur during Si shell burning and the 
early stages of collapse. However, as collapse ensues, 
quickly becomes large enough that /3-decay electrons are 
energetically blocked due to degeneracy. Similarly, j3'^- 
decay can also compete with electron-capture for nuclei 
with Qe+ > 2meC^, but in the core-collapse environment 
neutron-rich conditions are favored, and Qe+ is below 
this threshold. 

The importance of reactions mediated by the weak 
nuclear force, and specifically electron capture, as it per¬ 
tains to core-collapse was first demonstrated by Bethe 
et al. (1979). Not long after, the theory of stellar elec¬ 
tron capture was formalized by Fuller, Fowler, and New¬ 
man (FFN) (Fuller et al. 1982). In their pioneering work 
they published the first tabulation of weak interaction 
rates (/3^-decay and e^-capture) considering presuper¬ 
nova conditions where allowed Fermi and Gamow-Teller 
(GT) transitions dominate. Since then, advancements in 
computational resources have allowed for detailed nuclear 
shell-model calculations that have increased the accuracy 
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of the weak-interaction theory first outlined by FFN. Ma¬ 
jor weak-interaction rate tabulations that derive from a 
combination of experimental data and shell-model effec¬ 
tive interactions are the Oda et al. (1994) and Langanke 
& Martmez-Pinedo (2000) tabulations for sd- (A=17-39) 
and p/-shell (A=45-65) nuclei respectively. For heav¬ 
ier nuclei, the Shell Model Monte Carlo (SMMC) ap¬ 
proach has been employed to preserve nuclear properties 
such as the correlation energy scale in very large model 
spaces (Langanke et al. 2003). Langanke et al. (2003) 
have combined this method with an RPA technique to esti¬ 
mate electron-capture rates at densities and temperatures 
relevant during core-collapse for nuclei in the pfg/sdg- 
shell (A=65-112), which have come to be known as the 
LMSH rates. Juodagalvis et al. (2010) have also pro¬ 
duced a set of more than 2200 additional rates based 
on the same RPA technique but utilizing a Fermi-Dirac 
parameterization instead of the more computationally 
expensive SMMC calculations. The individual rates were 
not released, but instead these rates were averaged over 
NSE abundances and reported along a characteristic core¬ 
collapse {p,T,Ye) trajectory. Unfortunately, this is not 
suitable for the present study in which we investigate the 
detailed sensitivity on a nucleus by nucleus basis. 

In this work, we have implemented each of the rate 
tabulations listed in Table 1, which together contain 
445 rates for 304 unique nuclei over a large density 
and temperature grid. This library has been built as 
a standalone module and has also been implemented into 
the neutrino-interaction library NuLib (O’Connor 2015) 
for use in neutrino-transport routines employed by the 
spherically-symmetric, general-relativistic stellar collapse 
code GRID (O’Connor 2015)-see Section 4 for more infor¬ 
mation. Details on the density and temperature range 
for each of the included rate tabulations are shown in 
Table 1. The mass coverage of each rate table is shown 
in Figure 1. 

The LMP-I-LMSH rates were first implemented into 
a spherically-symmetric core-collapse simulation by Hix 
et al. (2003). They compared simulations with this set 
of shell-model based electron-capture rates against simu¬ 
lations that utilized the Bruenn (1985) prescription for 
electron capture. The evolution of the core-collapse phase 
and the structural differences in the core at bounce seen in 
that work were significant. In light of the differences that 
exist between theoretical estimates for electron-capture 
rates and those inferred from CE experiments, these re¬ 
sults motivate the need for a detailed sensitivity study. 

To handle the large number of nuclei not included in 
the tables, Hix et al. (2003) utilized an average electron- 
capture neutrino emissivity for all nuclei which lacked a 
shell-model based rate. Here, instead of performing aver¬ 
aging, we employ the approximate routine of Langanke 
et al. (2003), which is based on the parameterization of 
the electron-capture rate as a function of the ground state 
to ground state Q-value. This approximation was first 
described by Fuller et al. (1985) and was later parame¬ 
terized and fit to shell-model calculations in the p/-shell 
by Langanke et al. (2003). In this approximation, the 
electron-capture rate is written as: 


Aec = 


ln2 • B / T 


K 


[Fiiv) - “^xFsiv) + X^F2 {v)] 


( 1 ) 


and the neutrino-energy loss rate is, 

ln9 R / 'T \ ^ 

A.. = [F.ir,) - 2xF,{^) + x^F,{n)]^ 

( 2 ) 

where rUe is the electron mass, K = 6146 s, are Fermi 
integrals of rank k and degeneracy g, X = {Q ~ ^F)/T, 
g = y-f/ie/T, and T and pe are the temperature and elec¬ 
tron chemical potential. B (= 4.6) and AB (= 2.5 MeV) 
are fit parameters taken from Langanke et al. (2003) and 
respectively represent effective values for the transition 
strength and energy difference between final and initial 
excited states. 

In Figure 2 we compare the rate estimates from this 
approximation and those from the tables. As is easily 
seen from the figure, the variance of the shell-model rates 
depends sensitively on the density of the environment. At 
lower densities, where the electron chemical-potential and 
electron capture Q-value are comparable [pe ~ Qec), the 
location of excited states in the daughter nucleus, and the 
associated Gamow-Teller transition strength, sensitively 
determine the total electron-capture rate for a nucleus. 
New final states in the daughter nucleus become acces¬ 
sible as the electron Fermi-energy, which scales with 
increases beyond the energy required to populate them 
via allowed electron capture. The large scatter of the 
electron-capture rates at lower densities (Fig. 2) is be¬ 
cause the Fermi energy is comparable to the excited state 
energies of the daughter nuclei and the internal structure 
of each nucleus varies significantly. During collapse, as 
the density increases and the material becomes more neu¬ 
tron rich due to successive electron captures, both the 
magnitude of the average electron capture Q-value and pe 
increase. However, /Xg increases more quickly with density 
than the reaction Q-values do, and eventually pe S> Qec 
implying that the majority of the electron-capture chan¬ 
nels are open. In this regime the rate is less sensitive 
to the excitation energy spectrum of the daughter nu¬ 
cleus, and instead depends more strongly on the total GT 
strength across all possible final states. The decrease in 
the variance of the shell-model electron-capture rates in 
the higher density case of Figure 2 is a result of this. 

While the parameters of Eqs. 1 and 2 were originally fit 
from the LMP nuclei, there is reasonable agreement of the 
approximation with the other tabulated rates. Outside 
of these tables, significant deviations from the estimates 
of this approximation may exist, specifically for heavier 
neutron rich nuclei (Juodagalvis et al. 2008). But for the 
purpose of a sensitivity study, this approximation—from 
which the majority of the rates are calculated—is used 
as a base estimate off which the electron-capture rates 
may be varied. In Section 5.2 we show, given this set 
of rates, that heavier neutron-rich nuclei are the most 
important in the simulations. However, electron-capture 
rates developed from sophisticated theoretical models do 
not exist for individual nuclei in this region, and thus can¬ 
not be benchmarked against experimental measurements. 






4 


Sullivan et al. 


TABLE 1 

Density, temperature and mass ranges for the compiled weak rate set 


Model space 

Table 

s 

p 

pf pfg/sdg 

T (GK) 

LogiofpUe gem 3) 

Ref. 

FFN 

X 

X 

X 

0.01 - 100 

1.0 - 11 

Fuller et al. (1982) 

ODA 

X 

X 


0.01 - 30 

1.0 - 11 

Oda et al. (1994) 

LMP 

X 


X 

0.01 - 100 

1.0 - 11 

Langanke et al. (2003) 

LMSH 



X 

8.12 - 39.1 

9.22 - 12.4 

Hix et al. (2003); Langanke et al. (2001a) 

Approx. 

X 

X X 

X X 

- 

- 

Langanke et al. (2003) 


CD 

n 


c 

o 

o 

CL 
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Fig. 1. — Chart of the nuclear species included in each weak rate table. The table to which a species belongs is given by the color and 
legend in the figure. The Oda set contains rates for lower mass sd-shell nuclei (light blue), the LMP set contains rates for the intermediate 
mass p/-shell nuclei (green), and the LMSH set contains rates for the heavier mass pfg/sdg-shell nuclei near stability (red). The FFN 
tabulation provides rates across the sd and p/-shells (dark blue). Squares individually bordered in black are stable nuclei. The tables are 
mutually exclusive except for FFN which spans many nuclear shells. To distinguish between nuclei with rates from FFN and another table, 
the border of the FFN set has been outlined with a black and white line. 


Therefore, the estimates provided by the approximation 
of Eqs. 1 and 2 may be systematically off by a consider¬ 
able amount. As we will show, changes in the predicted 
rates for these nuclei have significant consequences for the 
simulations, motivating the need for experimental and 
theoretical efforts to constrain the rates of these species. 

3. CORE COLLAPSE AND THE ROLE OF ELECTRON 
CAPTURE 

Just prior to collapse, the temperature of the stellar core 
becomes high enough {T > 0.5 MeV) that the photon gas 
has sufficient energy to photodissociate nuclei into alpha 
particles and free nucleons. However, the density is also 
high (p > 10®gcm“^) resulting in large nuclear reaction 
rates that rapidly form nuclei from these light particles. 
The balance reached between these competing processes 
is known as Nuclear Statistical Equilibrium (NSE). If the 


entropy is sufficiently low and the mass fraction of free 
nucleons is small compared to that of nuclei, the most 
abundant nucleus in NSE is the species with the highest 
binding energy for a given electron-fraction, Ee(= Z/A) 
(Martfnez-Pinedo et al. 2000). A broad distribution of 
abundant nuclei forms due to finite temperatures which 
distribute the abundances around these peak nuclei. 

As collapse ensues and the central density increases 
through the first few decades, electron captures are the 
primary engine of deleptonization. Electrons are removed 
from the system and the produced electron-neutrinos 
(i^e) are able to freely stream out of the core, decreasing 
both Ye and total lepton fraction Y). As Ye decreases, 
peak abundances move toward neutron-rich nuclei, and 
the core begins to cool as Ve's carry away energy and 
entropy. Electron captures continue to dominate the 
neutrino transport during collapse until the last few mil- 
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Fig. 2.— Panel (a): Q-value dependence of electron-capture rates 
at two points along a core-collapse trajectory. The scattered points 
are tabulated (shell-model and SMMC) rates for each electron- 
capture reaction, while the black points are the approximate rates 
given by Eq. 1. Panel (b); The residual differences between logio of 
the shell-model rates and the approximate rates for each nucleus in 
the weak-rate library. An example residual is indicated on panel (a). 
When the density and temperature of a simulation evolve outside 
the range of the rate tables (see Table 1), rates are calculated via the 
approximate routines in order to avoid an artificial cut off imposed 
by the table boundaries. Rates are estimated between density and 
temperature grid points via monotonic cubic-spline interpolation 
as described by Steffen (1990). 

liseconds before core bounce. In these final moments, the 
central density reaches a few times 10^^ g cm“^, which 
is large enough that the neutrino mean free path begins 
to shorten due to coherent scattering on heavy nuclei. 
This increase in the i^e-scattering cross section results in 
a neutrino diffusion time that exceeds the collapse time, 
thereby trapping the electron neutrinos in the inward flow 
of matter. After this occurs, the conversion of electrons 
into electron-neutrinos via electron captures no longer 
removes leptons from the core. Instead, further electron 
captures increase the electron-neutrino fraction in 
order to conserve the now constant lepton fraction and 
bring the system of electrons and electron neutrinos into 
equilibrium. 

Prior to the work of Langanke et al. (2001a) it was 
believed that electron captures on free protons were of 
greater importance than captures on nuclei during col¬ 
lapse. The main considerations involved were that elec¬ 
tron capture on free protons has a higher rate owing to 
a smaller Q-value, and that nuclei with neutron number 
N > 40 have full p/-shell single particle states. Thus, 


the addition of another neutron via an allowed electron- 
capture transition would be Fermi-blocked. Langanke 
et al. (2001a) recognized that the many-body nuclear 
states have mixed configurations and do not follow a sim¬ 
ple Hartree-Fock filling of single particle orbitals. They 
also suggested that thermal excitation of nucleons to the 
( 7 g /2 orbital creates vacancies in the p/-shell, and together 
with configuration mixing, electron capture on bound 
protons is unblocked. Furthermore, because of the low 
entropy in the core, and the neutron-rich conditions, the 
abundance of heavy nuclei is several orders of magnitude 
higher than that of free protons, resulting in a higher 
overall electron-capture rate. Thus, because electron cap¬ 
tures on nuclei dominate, a detailed investigation of the 
contribution each species has to the deleptonization of 
the collapsing core, and in particular a determination of 
which nuclei are most important and therefore deserve 
further experimental and theoretical focus, would be of 
great value. 

4. CODES & METHODS 

4.1. NuLib 

In addition to the electron captures rates developed 
in this paper, other rates are needed to perform core¬ 
collapse simulations. The collection of rates we use is 
contained in NuLib (O’Connor 2015), an open-source, 
neutrino-interaction library available as a GitHub repos¬ 
itory at http://www.NuLib.org. NuLib contains rou¬ 
tines for calculating electron-type neutrino/antineutrino 
charged-current absorption opacities on nucleons with 
corrections for weak magnetism and nucleon recoil based 
on the formalism of Burrows et al. (2006) and Horowitz 
(2002). Neutrino emissivities for these processes are de¬ 
termined via Kirchhoff’s law which equates the absorp¬ 
tion rate of a equilibrium neutrino distribution to the 
emission rate of the underlying matter. Elastic scatter¬ 
ing of neutrinos on nucleons, and coherent scattering 
of neutrons on alpha particles and heavy nuclei is also 
included in NuLib. For the former we include correc¬ 
tions for weak magnetism and nucleon recoil, and for the 
latter we include corrections from ion-ion correlations 
(Horowitz 1997), electron polarization, and the nuclear 
form factor. Inelastic scattering of neutrinos on electrons 
is included based on the expressions of Bruenn (1985). 
Emissivities of heavy-lepton neutrino/antineutrino pairs 
via electron-positron annihilation and nucleon-nucleon 
Bremsstrahlung are computed ignoring final state neu¬ 
trino blocking. For neutrino-antineutrino annihilation, 
instead of computing the non-linear absorption opacity 
during the simulation, we make use of an effective absorp¬ 
tion opacity, which has been shown to be an excellent 
approximation for core-collapse supernovae (O’Connor 
2015). 

4.2. GRID 

We test our electron-capture rate implementation, and 
study the sensitivities of the core-collapse phase to 
these rates using the code GRID (O’Connor & Ott 2010; 
O’Connor 2015). GRID is an open-source spherically- 
symmetric general-relativistic neutrino-transport and hy¬ 
drodynamics code used for studying stellar collapse and 
the early stages of a core-collapse supernova. For details 
of the hydrodynamics module of GRID we refer the reader 
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to O’Connor & Ott (2010). The neutrino transport is 
handled though a general-relativistic, energy-dependent 
two-moment formalism for which extensive details can 
be found in O’Connor (2015). Our scheme numerically 
solves for the time evolution of the first two moments of 
the neutrino distribution function: the neutrino energy 
density and the neutrino momentum density. The sim¬ 
ulations utilize 18 energy groups logarithmically spaced 
between 0 and 250 MeV. Only electron type neutrinos 
are evolved until the central density reaches 10^^ gcm“^, 
after which electron anti-neutrinos and a characteristic 
heavy lepton neutrino are included. However, these lat¬ 
ter two neutrinos do not become important until core 
bounce has occurred. Spatial fluxes of the neutrino mo¬ 
ments are treated explicitly. Inelastic neutrino-electron 
scattering is handled explicitly until the central density 
reaches lO^^gcm”^ at which point an implicit treatment 
is used. Together NuLib and GRID provide a robust and 
extendable code base, making them ideal for the present 
study. 


4.3. Neutrino emission via electron capture 

Electron capture is associated with the emission of 
electron neutrinos and so the electron-capture rate is 
proportional to the integrated spectrum of Vg emitted per 
second. The rate for a particular nuclide, as tabulated in 
the implemented rate tables, is defined as the sum of the 
rates for each of the individual nuclear transitions 

A = (3) 

b 


where indices i and j correspond to levels in the par¬ 
ent and daughter nucleus respectively. The spectra of 
emitted neutrinos from the electrons capturing on nu¬ 
clei, described by the matter temperature T and electron 
chemical potential /ie, will vary based on the initial and 
final states involved owing to a different reaction Q-value, 

Qfp = Qg.s. +E,- E, (4) 

where Qg.s. is the atomic mass difference of the initial 
and final nuclei, and Ei and Ej are the excitation ener¬ 
gies of the populated states in the parent and daughter 
nucleus respectively. The most comprehensive solution 
to constructing neutrino spectra would be to coherently 
sum the spectra of neutrinos emitted from each nuclear 
transition. However this would rely upon rate tabulations 
for individual transitions which are not presently avail¬ 
able. Thus, we implement an effective neutrino spectra 
in terms of a single reaction Q-value, q, that is chosen to 
constrain the average energy of the spectrum to match 
that from the tabulated rates (Langanke et al. 2001b), 


n{E^,q) 


El{E,-qf 


N 


1 -I- expKE;, - q - pe)/kT} 


( 5 ) 


/ E^n{E^,q)dE, 

{EQ = ° 


Jn{E,,q)dE, Aec + A;3+ 


( 6 ) 



Pc (g cm-3) 


Fig. 3. — The average nuclear mass (divided by two), charge, 
and electron-capture rate versus central density for the three EOS 
utilized in this study. The colors indicate different EOS, while the 
line style indicate which quantity is plotted. All three EOS have 
nearly identical abundance distributions up to densities of g 

cm“'^. Beyond this point the TMA EOS has a heavier and slightly 
more neutron rich mass distribution compared to both SFHo and 
DD2, but maintains a comparable average electron capture rate 
overall. These simulations each utilize the sl5WW95 progenitor. 


where n{E^,q) is the neutrino distribution function and 
is normalized to the total electron-capture rate for a par¬ 
ticular nuclear species. Aj/, A^+, and Aec are the neutrino 
energy loss, positron emission, and electron-capture rates 
respectively. Eq. 6 is solved numerically for the effective 
Q-value, q = q^g, which then defines the effective neutrino 
spectrum for the electron-capture reaction of interest at 
a given p, T, and Yg. The approximate neutrino spectra 
generated in this way are unable to reproduce complex 
structure such as double peaking in the true neutrino 
distribution, which meLV occur when there is a resonant 
allowed transition {Qfp ~ 0) between an excited parent 
state and the daughter-nucleus ground state. However, it 
approximates singly-peaked neutrino distributions quite 
well (Langanke et al. 2001b). The spectrum is normalized 
to the total electron-capture rate via Gaussian-Legendre 
quadrature with an algorithm we developed to adaptively 
adjust the range of integration to the full width of the 
spectrum. 

Utilizing these spectra, the electron-capture neutrino 
emissivity for a given nuclear species is calculated as 

ciii.EQ — — Ej^n^ni^E}^^ q^g)^ (7) 

47r 

where the iZg’s are assumed to be emitted isotropically. 
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rij is the number density for the j-th nucleus, n{E}^^ q^g) 
is the neutrino spectra evaluated at the effective Q-value 
that solves Eq. 6 , and Ek indicates the energy of en¬ 
ergy group k. Evaluation of the emissivity is done point 
wise at the centroid of each energy bin, and has units 
erg/(cm^-sr-s-MeV). 

For estimates of the NSE number densities used above, 
we utilize several EOS from Hempel & Schaffner-Bielich 
(2010). In particular, we use the SFHo EOS and internally 
consistent NSE distribution developed by Steiner et al. 
(2013) for our primary EOS. We also compare with the 
DD2 (Typel et al. 2010) and TMA (Toki et al. 1995) EOS, 
each with self-consistent, but different, NSEs. The SFHo 
and DD2 EOS were chosen because they currently best 
satisfy both nuclear and astrophysical constraints (Fis¬ 
cher et al. 2014). Instead of meson self-interactions, the 
DD2 EOS implements density-dependent meson-nucleon 
couplings which have been used successfully to describe 
nuclear structure in a wide region of the nuclear chart 
and have also been tested in heavy-ion collisions (Typel 
et al. 2010). For nuclear masses, the SFHo and DD2 
NSE distributions rely on the Finite Range Drop Model 
(FRDM) from Moller et al. (1995) and Moller et al. (1997), 
whereas the TMA EOS utilizes a mass table calculated 
by Geng et al. (2005). For consistency, in addition to NSE 
abundances, these mass distributions are also utilized in 
the calculation of reaction Q-values for use in Eq. 1. 

5. SENSITIVITY STUDY 
5.1. Reference simulations 

In order to establish reference simulations off which 
variations are performed, we have utilized the widely 
studied 15 solar mass, solar metallicity progenitor 
sl5WW95 (Woosley & Weaver 1995), as well as sl2, s20, 
and s40 from WH07 (Woosley & Heger 2007), which span 
the range of stellar compactness ^ 2 . 5 (O’Connor & Ott 
2011) in this model set. For more details, see the progen¬ 
itor sensitivity subsection below. For each progenitor we 
use the SFHo EOS described above, and in addition, for 
S15WW95 we also employ the DD2 and TMA EOS and 
NSE distributions. 

For each reference simulation, a full complement of 
neutrino-interaction microphysics is incorporated via 
NuLib, which includes the newly implemented weak rates 
library described here. The weak-rate tables were in¬ 
cluded using the following priority hierarchy: LMP > 
LMSEl > Oda > Approx., ensuring that rates from 
sources with higher priority are utilized where rate es¬ 
timates from multiple sources exist. Approx, indicates 
the parameterized rate approximation of Eqs. 1 and 2, 
which is used for nuclei not included in the tables and for 
regions of density and temperature which are beyond the 
limits found in Table 1. For consistency, only tables that 
derive from shell-model calculations are utilized. 

For each progenitor and EOS we perform collapse sim¬ 
ulations in GRID and follow the evolution until at least 
^ 100 ms after bounce. The collapse proceeds as described 
in Section 3. Differences in the collapse evolution for dif¬ 
ferent progenitors stem from the hydrostatic conditions 
in the cores of these massive stars at the onset of collapse. 
For stars with large ^ 2 . 5 ; larger central temperatures are 
needed to balance gravity. This gives lower central densi¬ 
ties, and therefore less electron capture during the final 


stages of stellar evolution. The range of initial central Yg 
goes from ~0.422 for the sl5WW95 model to ^0.447 for 
the s40WH07 model, or a range of ~ 6 %. After neutrino 
trapping sets in, we see a range of trapped lepton-fraction 
of ~0.288-0.297, where s40WH07 and sl2WH07 have the 
minimum and maximum trapped Yj, respectively. The 
overall higher deleptonization rate for the more compact 
progenitors is due to both longer collapse times and larger 
matter temperatures, which enhance the electron-capture 
rates. 

Simulations utilizing different EOS, while holding all 
else constant, demonstrate only small variations in the 
density, temperature, and Ye central-zone trajectories up 
to bounce. Figure 3 details the abundance distributions 
for each EOS, as well as the resulting average electron 
capture rate along a collapse trajectory. The NSE distri¬ 
butions of all three EOS are largely similar early on, but 
differences in the mass table of the TMA EOS cause it to 
diverge from the others starting around 1 — 2 x lO^^gcm”^. 
However, differences are seen in the electron-capture rate 
only after central densities of 2 x lO^^gcm”^, where any 
effect on the evolution is suppressed because of neutrino 
trapping. Near nuclear saturation density, however, the 
differences in EOS begin to play a more important role. 
The density-dependent couplings of the 01)2 EOS, for 
instance, result in higher central temperatures at bounce. 
However, since the average rate in simulations utilizing 
each of the EOS are nearly identical, they result in a 
difference of trapped lepton-fraction of only a fraction of 
a percent. For more information on the sensitivity of the 
electron and lepton fractions to the EOS during collapse 
see Fischer et al. (2014). 

Together, these reference calculations span a wide 
range of progenitor and EOS dependences that ensure a 
configuration-independent assessment of the core-collapse 
sensitivity to electron capture on nuclei, and furthermore 
demonstrate the universality of collapse. In what fol¬ 
lows, we detail the results primarily for variations on the 
sl5WW95+SFHo reference simulation, but describe any 
significant differences, where they exist, in relation to 
variations on the other progenitor-l-EOS reference simula¬ 
tions. 

5.2. Species dependent sensitivity 

To understand the sensitivity of core-collapse to differ¬ 
ent regions of electron capturing nuclei, we use the central 
zone collapse profile from our reference simulation to de¬ 
compose the change of the electron fraction with time, Yg, 
into the electron captures of each nuclear species. While 
using only the central zone is an approximation, we jus¬ 
tify this by noting the observation by Liebendorfer (2005) 
that the electron fraction profiles typically correlate quite 
well with density during the collapse phase. Therefore, 
matter will generally have the same electron capture his¬ 
tory. We calculate Yg and account for iZg re-absorption in 
the following energy-dependent way, 

• i dTTQ! ^ Acfe • ?7,(efe) ( Ek\ 

‘‘W.E —S- 

where Yf is the time derivative of the electron fraction due 
to electron captures on the Ah nuclear species, a accounts 
for the general relativistic time dilation, Na is Avogadro’s 
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constant, p is the density, is the energy of the fcth energy 
bin, Ae/j is the kth energy bin width, rji is the emissivity 
of species i and 1 — ^ is the neutrino blocking factor that 
accounts for re-absorption as collapse approaches weak 
equilibrium. Along with a hydrodynamical correction 
due to advection of electrons into the central zone, the 
time integral of Eq. 8 added for all nuclei reproduces the 
full time dependent Ye profile of the central zone during 
collapse, indicating that electron captures on heavy nuclei 
singularly drive the deleptonization of the central zone. 

With this method, the deleptonization history due to 
each nucleus can be individually investigated. At these 
densities and temperatures, NSE diversifies the abundant 
nuclei, ensuring that no single nucleus dominates the 
deleptonization. There are, however, subsets of nuclei 
that contribute more than others to the reduction of Yg. 
A nuclear-mass dependence can be studied by binning the 
the contribution to \Yf.\ from each nuclide into nuclear 
mass bins and tracking the evolution of each region up 
to neutrino trapping. Figure 4a plots the deleptonization 
rate in the core for different nuclear mass bins, as the 
central Yf, progresses from its progenitor value to its value 
when weak equilibrium is achieved, just prior to bounce. 
Early on, before the collapse becomes strongly dynamical, 
nuclei in both the mass range 25<A<65 {sd+pf-shell) 
and those in the 65<A<105 (pfg/sdg-shell) comprise the 
main component of the deleptonization. However, during 
the strongest push toward neutron-rich conditions, where 
Ye rapidly changes from ^0.41 to ^0.28, nuclei with mass 
A >65 dominate the evolution as seen by the red and light 
blue curves in Figure 4a. Unfortunately, the most precise 
electron-capture rate estimates fall below this region and 
instead, the rates are set primarily by the approximation 
of Eq. 1. 

It is also useful to understand the specific nuclei that 
have the largest integrated contribution to core delep¬ 
tonization up to neutrino trapping. In panel (b) of Fig¬ 
ure 4 we plot the 500 nuclei with the largest integrated 
I Ye I from t = 0 to the trapping time—when densities are 
in excess of 2 • 10^^ g cm“^. This reveals the channel 
through which the bulk of electron captures operate. The 
central electron-fraction at the trapping density is repro¬ 
duced by subtracting the sum of this quantity over all 
nuclear species from the initial electron-fraction, 

Ye{t = ^trapping) - Ye{t = 0) - ^ AY^, (9) 

i 

where AY^ is shown in Figure 4b, and the component 
of Ye due to advection of electrons into the central zone 
(otherwise making this relation exact) is left out for sim¬ 
plicity. Within the pfg and sdg-shells we find that pri¬ 
mary contributors to the deleptonization phase of collapse 
are neutron rich nuclei near the N=50 and N=82 closed 
neutron shells. 

To confirm these results, we also gauge the sensitiv¬ 
ity of the collapse phase to localized groups of nuclei 
by employing a statistical resampling technique where 
sets of nuclei are removed from the simulation. This 
method is based on well known statistical resampling 
methods such as bootstrap and jackknife resampling (Wu 
1986). Specifically, a rectangular region centered on a 
nucleus and spanning all nuclei within ±3 isobars and 


t-tb (ms) 

<-100 - 50 - 40-30 -20 -10 



Fig. 4.— (a) The contribution of nuclear electron capture to the 
change of the matter electron-fraction with time. The contours are 
the binned sums of l^el for each species in several mass regions. For 
reference, the central density at t = -20, -10, -5, -2, and -1 ms 
is 1.41 ■ 10i\ 4.06 ■ 10^^ 1.42 ■ 10^^, 8.49 ■ lO^^, 3.30 ■ 10^^ g cm-^ 
respectively, (b) Top 500 electron capturing nuclei with the largest 
absolute change to the electron fraction up to neutrino trapping. 
The color scale indicates \Ye\ integrated up to the trapping time, 
occurring when pc 2 • 10^^ g cm~^, such that the total electron- 
fraction at this point is equal to its initial value less the sum of AVe, 
the plotted quantity, over all nuclides. Calculations are based on 
the sl5WW95-hSFHo reference simulation. The rectangular outline 
indicates the size of the sampling region used in the statistical 
resampling study, and also the set of nuclei which exhibited the 
largest changes to the simulations when excluded from the electron- 
capture calculations. 

zb5 isobaric chains is removed from the calculation of the 
electron-capture neutrino emissivity. An example of such 
a removed region is drawn on Figure 4b. This technique 
is employed in 48 simulations with resampling performed 
uniformly across the nuclear chart. Using this technique 
we find the simulations are most sensitive to nuclei in 
the mass range 74-84 with Z/A (= Yg) between 0.36-0.44, 
corresponding to nuclei near ^^Ni, ^^Cu, and ^^Zn. These 
results agree with the calculations performed above, 
and indicate that species near the N=50 magic number 
have the largest contribution in magnitude to the change 
in the electron fraction overall. The impact of remov¬ 
ing these species from the simulation corresponded to a 
change of inner-core mass at bounce of > 10%, whereas 
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Fig. 5. — Projection of the (trapped) lepton fraction at = 
4 ■ 10^^ gcm“® as a function of the electron-capture rate scaling 
factor for progenitor-|-EOS reference simulation. In all the cases 
the lepton fraction begins to increase if the capture rate becomes 
too high because of a dramatic increase in the electron neutrino 
absorption cross section. The asymmetry seen here indicates that 
those quantities which depend on Yij^ are likely to be more sensitive 
to a reduction of the electron-capture rates due to a systematic 
overestimate in the base rates, than they are to an increase due to 
an underestimate. 

resampling in other regions resulted in variations of only 
a few percent. 

The electron-capture rates for these nuclides rely en¬ 
tirely on the approximation of Eqs. 1 and 2, which were 
fit originally to rates of lower-mass mid-shell nuclei near 
stability. Therefore, in the region indicated by the above 
two studies, the approximation is largely uncertain and 
may be systematically off by a significant amount. For 
instance, these estimates do not account for nuclear struc¬ 
ture effects that may occur near the N=50 closed neutron 
shell. Depending on the nuclear configurations, ther¬ 
mal excitations, and increasing dependence on forbidden 
transitions, Pauli blocking may considerably reduce the 
electron-capture rates in this area. Given that the change 
of inner-core mass at bounce was largest when the rates of 
these nuclides were decreased to zero as compared to any 
other set, and that without any evaluative measurements 
the uncertainties in these rates remain large, experimental 
and theoretical work should focus here. Any substantial 
changes to the electron-capture rate estimates for these 
nuclei will likely have a relatively large impact on simula¬ 
tion predictions for the PNS formation, and will therefore 
help to constrain important collapse and pre-explosion 
phase quantities. 

5.3. Systematic variations 

To study the strongest impact of variations in the 
electron-capture rates, we perform simulations in which 
the rate for each A>4 nuclide is systematically scaled by 
factors of 10, 4, 2, 0.5, 0.25, and 0.1. In this way, the 
structure of the rates as seen in Figure 2 is preserved 
(the lower panel of residuals is unaffected), but the dis- 
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Fig. 6. — Comparison of the central electron, electron-neutrino 
(left) and lepton (right) fractions in which the nuclear electron- 
capture rate for every species has been scaled by factors shown in 
the legend. Warmer colors indicate a higher overall electron capture 
rate, and cooler colors indicate a lower rate. The dashed black line 
indicates the reference sl5WW95-|-SFHo simulation. 


tribution of rates is shifted to larger or smaller values 
depending on the scaling factor. Systematic shifts of the 
rates emphasize the role of electron capture as a regulator 
for entropy and temperature in the simulations. By in¬ 
creasing the rates, more neutrinos are emitted and escape 
during the initial stages of collapse, thereby increasing the 
evaporative neutrino-cooling. Furthermore, because the 
dominant source of matter pressure is electron degeneracy, 
increased electron-capture rates accelerate the collapse. 
This impacts the matter profiles outside the shock in the 
early post-bounce phase. Decreasing the rates has the 
opposite effect, the entropy, temperature, and electron 
fraction of the core are significantly higher because less 
cooling takes place. 

The evolution prior to and right at pc = 2 • 10^^ gcm“^ 
(which is the density that defines neutrino trapping) is 
what sets the final value of the trapped lepton and elec¬ 
tron fractions, which are important due to their direct 
impact on the formation of the PNS. For all the reference 
simulations we observe a minimum in the trapped lepton 
fraction occurring for a systematic scaling factor of ap¬ 
proximately four. The minimum that forms can be seen 
in Figure 5. Scaling by ten slightly reverses the downward 
trend, and increases the trapped lepton-fraction from its 
minimum value. This behavior is the result of electron- 
neutrino capture on heavy nuclei becoming the primary 
source of opacity, exceeding what is typical as a result 
of coherent i/g-scattering. When the rates have been en¬ 
hanced by a factor of ten, the ratio of the absorption and 
scattering opacities, Ka/^'si surpasses unity already by 
central densities of 3 • 10^^ g cm“^ and Kq ~ 4ks by the 
time Pc = lO^^g cm“^. Absorption cross sections are then 
large enough to trigger an early onset of neutrino trapping 
at densities lower than what is found for the reference 
rates. The consequence is that electron capture has a 
smaller window of deleptonization, leading ultimately to 
a higher Yj overall. 

The range of electron fractions near core bounce is com¬ 
mensurate with the range of trapped lepton-fractions so 
far described, see Figure 6. As mentioned above, varia¬ 
tions of Ye (and Yi) on this level are of importance due to 
its direct impact on the formation of the PNS and the su- 
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Fig. 7. — The electron fraction, entropy, density, and velocity as a function of enclosed mass at four times during a core-collapse simulation, 
spanning 6 ms around bounce, including the collapse phase just after the onset of neutrino trapping (a), core bounce (b), and 1 ms (c) and 5 
ms (d) after bounce during which the shock has begun its outward trajectory. The reference simulation (sl5WW95-|-SFHo) is shown in 
dashed black. Simulations shown in color have electron-capture rates scaled systematically for all species by factors of 10, 4, 2, 0.5, 0.25, and 
0.1. Core bounce {t — ti, = 0) shown in panel (b) is defined as when the entropy at the shock front exceeds 3.0 kj/baryon. 


pernova shock. Electron fraction, entropy, density and ve¬ 
locity profiles are shown in Figure 7 for sl5WW95-|-SFHo 
at -1, 0, 1, and 5 milliseconds relative to bounce. Of par¬ 
ticular interest, we find that the mass of the forming PNS 
inner-core at bounce, seen as the mass behind the steep 
velocity gradient in panel (b), varies on the order of ^0.1 
Mq, and up to ~0.2 PIq five milliseconds after bounce. 
The asymmetry observed in the trapped-lepton fraction, 
where scaling the rates by 0.1 had a more dramatic effect 
than scaling by 10, translates directly to the variation 
of the inner-core mass at bounce (-|-16/-4% from the 
reference). The result we find is that the forming PNS 
has a lower bound on the inner-core mass at bounce over 
the range of electron-capture rates explored. Because 
the rates are already high, and therefore the absorption 
opacity is already almost comparable to the scattering 
opacity, the range of inner-core mass at bounce comes 
mainly from simulations with decreased rates relative to 
the base simulation. 

In addition to the direct impact on core dynamics and 


structure, the neutrino emission at bounce is found to 
be very sensitive to these variations. Figure 8 shows 
the neutrino luminosity 500 km from the center for the 
different neutrino species as a function of time. Prior 
to bounce the i/e-luminosity begins to rise from electron 
captures on bound protons in nuclei, but is quickly regu¬ 
lated by neutrino trapping, causing a down turn in the 
luminosity. During this time the core is very sensitivite 
to the nuclear electron-capture rates as the entropy is low 
enough that heavy nuclei dominate the available mass. 
Scaling the rates for each nucleus using the same system¬ 
atic factors results in a 40% variation of the i^g-luminosity 
before bounce. During bounce, the electron-neutrino 
burst—seen as the peak luminosity in the left panel of 
Figure 8 —is powered primarily by electron capture on free 
protons. The core-bounce and shock liberates nucleons 
from their bound states and the entropy rises causing a 
signihcant increase in the nucleon and light particle abun¬ 
dances. That said, while the electron-capture rate on free 
protons, is not adjusted in these simulations, a range 
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Fig. 8. — The neutrino luminosity as measured at a radius of 
500 km as a function of time after bounce in our sl5WW95-|-SFHo 
simulation set. Electron-capture rate scaling factors are shown in 
the legend, where contours with warmer colors have higher rates, and 
cooler colors have lower rates. While the peak electron-neutrino 
luminosity is considered particularly stable across core-collapse 
simulations, it varies significantly with variations of the electron- 
capture rates on medium-heavy nuclei. When the rates are at 
their lowest (xO.l case), the shock reaches the neutrinosphere more 
quickly than in the other simulations. This results in a larger 
luminosity in the peak electron-neutrino burst because more i^eS 
are able stream out of the core at early times. The opposite is true 
when the rates are higher, the neutrinosphere and shock converge 
much more slowly, and so the neutrinos spend more time diffusing 
out of the inner core, reducing the peak luminosity but distributing 
it out to later times. 


Sensitivity to 

Progenitor variations Eiectron capture rate variations 



of ±20% relative to the reference peak t/e-luminosity is 
observed. 

We find that these dramatic variations of the peak 
electron neutrino luminosity are a result of alterations 
to the neutrinosphere and shock convergence-timescale. 
Specifically, when electron captures on nuclei are weaker 
(scaling by 0.1), the inner-core mass that forms at bounce 
is significantly larger. This results in more kinetic en¬ 
ergy transferred to the shock, allowing it to sweep up 
mass more quickly. In Figure 7 this can be seen by the 
broadening of the distribution of shock locations in mass 
between the different simulations in the velocity plot 5 
ms after bounce (bottom-right) as compared to t — tb = 
0. Also, with a weaker overall rate the opacity will be 
lower, allowing the neutrinosphere to move in to lower 
radii more quickly. The combination of these effects result 
in the shock and neutrinosphere radii converging earlier 
for the simulations with lower electron-capture rates, and 
later for simulations with higher rates, up to a difference 
on the order of 3.5 ms. Thus, electron capture on protons 
liberated by the shock produce neutrinos that are able to 
reach the neutrinosphere earlier and freely stream away, 
contributing to a larger peak luminosity when the nu¬ 
clear electron capture rate is systematically lower. On 
the other hand, when the nuclear electron capture rate is 
high, the emitted neutrinos diffuse more slowly through 
the core, and reach the neutrinosphere at later times, thus 
strongly quenching the peak luminosity but spreading out 
the emission to later times. Due to the high luminosity 
of the electron-neutrino burst near the time of bounce, it 


Fig. 9. — The full range of sensitivity of the PNS inner-core mass, 
central entropy, and central temperature at bounce as well as the 
peak i^e-luminosity, the peak average Ue energy, and the average 
i/e energy prior to neutrino trapping, owing to variations of the 
progenitor model and electron-capture rates. Thirty two progenitors 
were utilized from the WH07 model set of Woosley Sz Heger (2007) 
for producing the progenitor bars (red) in the figure. Each bar of 
the electron-capture rate variations derives from simulations where 
the rates have been systematically scaled by factors of 10, 4, 2, 
0.5, 0.25, and 0.1. The horizontal tick represents the value of the 
reference simulation for the tested Progenitor -h EOS combination. 
The window ranges are chosen so that the progenitor sensitivity 
bars are of equal size across each of the plotted parameters. 

is a candidate for detection from a galactic core-collapse 
supernovae in Earth-based detectors sensitive to electron 
neutrinos, e.g. those with a detector volume composed 
of liquid Argon. And while such measurements are not 
presently of high enough precision to resolve each vari¬ 
ation seen here, they may indicate the total amount of 
electron capture occurring at core bounce. 

5 . 3 . 1 . Progenitor model sensitivity 

In order to evaluate the significance of the electron- 
capture systematic sensitivity studies described so far, we 
have chosen to test against a study of the progenitor depen¬ 
dence of the core-collapse phase. Drawing from the larger 
set of progenitors from which the reference progenitors of 
the electron-capture study belong, the 2007 non-rotating 
solar-metallicity single-star model set from the stellar 
evolution code KEPLER (Woosley & Heger 2007) was 
utilized. This model set contains the presupernova con¬ 
figuration of 32 stars ranging in zero-age-main-sequence 
(ZAMS) mass from 12 Mq to 120 Mq- s 12WH07 and 
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sl20WH07 respectively. Simulations of these progenitors 
exhibit a ~3.5% range of trapped lepton-fraction (0.288 - 
0.298), a ^4% range of inner-core mass at bounce (0.473- 
O. 49 IM 0 ), and a ^9% range of electron-neutrino peak 
luminosity (5.19 — 5.65 • 10®^ erg s“^) during the neutrino 
flash occurring just after core bounce. 

Figure 9 compares the progenitor model and electron- 
capture rate dependence of several structural and neutrino 
quantities during collapse. The range of inner-core mass 
and peak v^, luminosity seen from employing the WH07 
progenitor model set are each approximately a factor of 
5 smaller than the ranges seen from varying the electron- 
capture rates across all progenitor-|-EOS references. On 
the other hand, the range of central entropies and tem¬ 
peratures at bounce are comparable between the two 
sensitivity studies. The v^. average energies just prior to 
neutrino trapping and during the deleptonization burst 
are also compared in Fig. 9. The neutrinos emitted during 
the luminous burst just following core bounce are of higher 
energy than those emitted earlier because they arise pri¬ 
marily from electron capture on free protons. They also 
decouple from the core at a much hotter and denser neu- 
trinosphere than prior to bounce, yielding higher energy 
neutrinos. In both of the sensitivity studies, variations of 
the electron-capture rates and of the initial stellar models, 
we find that the range of average neutrino energy during 
peak emission is comparable (« ±0.5 MeV). 

While captures on free protons contribute only 
marginally to deleptonization in the central zone, fur¬ 
ther out in the iron core, where the densities are lower 
(and Ye’s are higher), electron captures on protons con¬ 
tribute to the deleptonization, especially in cases where 
we suppress electron captures on nuclei. The capture 
of electrons on these free protons produces neutrinos of 
a higher average energy, commensurate with the large 
spread seen in the bottom panel of Fig. 9 (which is taken 
when the central density is 3 x 10^^ gcm“^, but present 
from the onset of collapse). Another contribution to the 
energy spread is the systematic shift of electron captures 
to more neutron rich nuclei as the electron-capture rates 
are increased and the matter becomes more neutron rich. 
These neutron-rich nuclei have more negative Q-values, 
yielding lower energy neutrino emission. Both of these 
effects result in a dispersion of average neutrino energies 
early on that is several factors larger than what is seen 
in the progenitor simulations. 

Finally, we note that while the peak luminosity is only 
weakly dependent on the progenitor model, the post¬ 
bounce pre-explosion luminosity of all six neutrino species 
have strong progenitor dependences (O’Connor & Ott 
2013), and we find these pre-explosion luminosities are 
much less sensitive to the nuclear electron-capture rates 
comparatively-see panel (b) of Figure 8. The diverging 
of the luminosities seen at t — tb = 120 ms is due only to 
the difference in collapse times between the simulations 
which carries over to the evolution of the mass accretion 
rate after bounce. 

5.4. Monte-Carlo variations 

In addition to the possibility of systematic errors in the 
electron-capture rates, we also explore the effect of statis¬ 
tically distributed variations. Such an investigation is of 
great importance if the effect of an approximation such as 


Eq. 1 is to be understood. The main flaw in a continuous 
function for rate estimation across many nuclear species 
is the loss of structure, which would otherwise serve to 
statistically distribute the rates on a reaction by reaction 
basis (see Fig 2). To study this effect we performed a 
Monte-Carlo (MC) variation of the electron-capture rates. 
Using an analytic description of the electron-capture rate 
distribtutions, such as a Gaussian or Poisson distribu¬ 
tion, is likely to be inaccurate. Instead, we MC adjust 
the approximate rate for each species by adding to its 
logio(AEc) a value randomly chosen from a distribution 
created from the residuals of the tabulated rates and the 
approximate rates, i.e. 

logio(AEc“®) - logio(AEc'^'^) (10) 

where i is an index running over all the tabulated reac¬ 
tions. 

In constructing this distribution, it is important also 
to preserve the Q-value dependence of the residuals that 
can be seen in Figure 2b. We do this by seperating the 
residual distribution into subsets so that the reaction-rate 
residuals in each subset have similar Q-values. To do 
so we have chosen a Q-value binning of 2.5 MeV, but 
have also tested binnings of 5.0 MeV and 10.0 MeV which 
resolve the Q-dependence less, but have more counts per 
bin from which to sample. With this method we MC 
generate pseudo electron-capture rates that retain the 
Q-value dependence of Eq. 1, but statistically distribute 
the approximate rate according to the variance of the 
rates calculated in the shell-model. Seven simulations for 
each binning were performed. 

As mentioned before, at low densities the electron- 
capture rate depends strongly on the energy levels of 
the initial and final nuclides because the electron chemi¬ 
cal potential is comparable to the excitation energies of 
the allowed Gamow-Teller transitions. As the electron 
chemical potential increases, it encompasses a larger range 
of excitation energies which results in the electron-capture 
rate becoming sensitive primarily to the total strength. 
In the low-density case of Figure 2a the approximation 
of Eq. 1 while appearing to decently reproduce the mean 
of the shell-model rates, actually has a mean approxi¬ 
mately a factor of two lower than the tabulated rates. As 
the density increases, this difference between the mean 
electron-capture rate estimated by the approximation and 
the shell-model rates decrease. Thus, the approximation 
better reproduces the mean rate in the high density case 
of Figure 2a. Because our MC simulations are based 
on residual distributions of the tabulated rates and the 
approximation, the average rate produced in each MC 
trial also has this bias. 

In Figure 10 we plot the min-to-max band representing 
the range of lepton fraction observed from all of the MC 
simulations. The band drawn corresponds to a 2.5 MeV 
binning of the residual distributions from which the MC 
sampling was performed. For the reasons just described, 
the band has lower electron and lepton fractions than the 
reference at low densities, but becomes more statistically 
distributed around the reference at higher densities, near 
5T0^^g cm“^. The lepton fraction band width varies 
from about a half percent initially, to its largest value of 
^2.5% just before neutrino trapping, and then decreases 
back to ~1.5% before bounce. Altogether, we observe no 
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Fig. 10. — Lepton fraction band as a function of central den¬ 
sity for the Monte Carlo (MC) described in the text. The band 
represents data points from seven MC simulations based on the 
sl5WW95-|-SFHo reference, where the band width represents the 
min to max range of Yj. The rates are varied by drawing from 
the Q-value dependent residual function shown in Figure 2 after 
binning it into 2.5 MeV energy bins. This binning was chosen as it 
best tracks the Q-value dependence of the residual distribution. 

significant impact on the core dynamics or the neutrino 
transport and therefore conclude that any statistically 
distributed scatter in the estimations of the electron- 
capture rates, such as those seen in Figure 2, will likely 
not impact the models. 

6. CONCLUSION 

Nuclear electron capture has long been understood to 
play an important role in the dynamics of core-collapse 
supernovae and large efforts have been undertaken to 
produce reliable estimates of electron-capture rates for 
astrophysical contexts. Although significant progress has 
been made in benchmarking theoretical electron-capture 
rates by comparison with charge-exchange experiments 
(especially using shell-model calculations) (Cole et al. 
2012 ), large uncertainties remain for neutron-rich nuclei 
and nuclei beyond A=65. Furthermore, sophisticated 
shell-model estimates for electron-capture rates exist only 
for a small subset of the large number of nuclei that con¬ 
tribute strongly. In this work we explore the implications 
of uncertainties in the electron-capture rate estimates for 
the core-collapse and early post-bounce phases of fully self 
consistent, general relativistic, core-collapse supernova 
simulations with comprehensive neutrino transport. 

6.1. Most important nuclei 

For the reference simulation we calculate the contri¬ 
bution of each nucleus to core deleptonization and also 


perform a statistical resampling study, both of which 
identify those species whose rate should be known most 
precisely due to their significance in the simulations. With 
the given set of electron capture rates—from shell-model 
estimates to the approximate estimates of Eqs. 1 and 2— 
we find that the simulations are most sensitive to neutron 
rich nuclei in the upper pf and pfg/sdg-sheWs. 

Specifically, in these simulations nuclei near the A~80, 
N^50 closed neutron shell contribute the bulk of core 
deleptonization, and when removed from the simulations 
result in noticeable changes to the protoneutron star for¬ 
mation, with a significantly larger impact than when any 
other group of nuclei are removed. However, because so¬ 
phisticated estimates from nuclear theory are not available 
for individual nuclei in this region, the electron-capture 
rates for these species have been accounted for in the 
past via simple averaging techniques and in this work via 
an approximation that has been fit to the LMP rate set. 
While this approximation reasonably reproduces the aver¬ 
age electron capture rate for sd and pf shell nuclei near 
stability, rates for heavier neutron-rich nuclides will likely 
diverge from what is predicted by this parameterization. 

6 .2. Impact of uncertainties 

We evaluate the impact such uncertainties may have 
by varying the electron-capture rates for more than 6000 
nuclei statistically, about the approximate prediction, and 
also systematically. On one hand, we find that statisti¬ 
cal variations of electron-capture rates effect the overall 
dynamics and neutrino emission only weakly, producing 
marginal changes to the simulations. These findings indi¬ 
cate that the lack of structural variation that distributes 
the rate estimates from one species to the next is not 
crucial to the simulations. 

On the other hand, the average electron capture rate 
across a region of nuclei strongly determines the overall 
impact of those constituent nuclei. By systematically 
varying the electron-capture rates by factors between 10 
and 0.1 we observe dramatic variations in the inner-core 
mass (-|-16/-4%) and the electron-neutrino luminosity (± 
20%) at and near bounce, respectively. Comparing with 
32 simulations utilizing different progenitor models, we 
find that this range of inner-core mass and peak neutrino- 
luminosity is 5 times as large as that seen when varying 
the progenitor models. 

We also find that the nuclear electron-capture rates are 
already large enough in the reference simulations that in¬ 
creasing them beyond their base values has a considerably 
smaller effect than decreasing them. This has compelling 
implications. Rates for A^80 nuclei near the N=50 shell 
gap, which have been shown in this work to be the primary 
contributors to the overall impact of electron captures 
during core collapse, may be overestimated by Eq. 1 due 
to Eermi-blocking at the closed neutron shell. Combined 
with a greater overall sensitivity to the systematic de¬ 
crease in electron-capture rates, changes to the collapse 
and early post-bounce phases of the simulations may be 
as significant as those seen in this study if the current 
rates of these nuclei are found to be overestimated. 

6.3. Goals for future studies 

For these reasons, it is important that experimental 
and theoretical efforts be aimed at nuclei which span 
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the region on the chart of isotopes between stability and 
the neutron drip-line in both the pf and pfg/sdg model 
spaces, and further expand on the work that has been 
carried out for (near-)stable nuclei in the p/-shell. Since 
data from (n,p)-type charge-exchange experiments for 
nuclei in the pfg/sdg-sheW and for neutron-rich nuclei 
in the pf and pfg/sdg-sheW are scarce, new experiments 
are required to obtain a sufficient set of data to bench¬ 
mark current and future theoretical estimates. To this 
end, presently feasible experiments on neutron-rich nu¬ 
clei at and near stability with 60<A<120 should add to 
the few cases that have been measured in this region. 
With the higher beam intensities that will be available at 
next generation rare isotope facilities, future experimental 
programs should focus on the neutron-rich component 
of the primary electron-capture channel shown in Fig¬ 
ure 4b. In particular, investigation of nuclei in the A^80 
and N~50 region should take precedence, as changes to 
their electron-capture rates will significantly constrain 
the core-collapse dependence on nuclear electron-capture. 
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